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1 .  Introduction 

This  final  technical  repon  concerns  research  performed  over  a  three-year  period 
ending  February  14, 1997. 

Researchers  who  participated  in  the  work,  in  addition  to  the  author,  were  Dr.  Jong- 
Soo  Kim,  Professor  Paul  Clavin  and  Mr.  Chae-Hoon  Sohn.  Professor  Clavin  was 
concerned,  in  particular,  with  turbulence  and  noise  effects  on  combustion  instabilities.  Mr. 
Sohn  was  concerned,  in  particular,  with  droplet  combustion  as  an  amplification 
mechanism,  especially  LOX  droplets  in  gaseous  hydrogen  atmospheres.  Dr.  Kim  and  the 
author  were  involved  in  all  aspects  of  the  work,  including  those  of  Professor  Clavin  and 
Mr.  Sohn.  The  following  sections  summarize  all  of  this  research. 

2 .  Publications 

All  of  this  work  is  now  published  in  the  open  scientific  literature.  References  [1] 
through  [6]  list  these  publications.  The  reader  is  referred  to  these  publications  for  details  of 
the  research.  The  following  sections  therefore  merely  provide  summaries  of  the  research. 

3 .  Theoretical  Explanation  of  Acoustic  Instability  Boundaries 

It  has  been  known  for  some  time  that  the  most  severe  instabilities  in  liquid  rockets 
are  acoustic,  that  is.  sound  waves  may  be  amplified  by  interactions  with  combustion 
processes  to  an  intensity  high  enough  to  damage  the  engine.  However,  the  acoustic 
instability  mechanisms  have  been  poorly  understood,  so  that  as  a  design  guideline 
engineers  had  to  rely  on  empirical  methodologies  based  on  previous  engine-test  results  to 
estimate  where  the  boundaries  of  acoustic  instability  lie.  Good  understanding  of  the 
mechanisms  by  which  such  empirical  relations  emerge  could  help  in  the  design  of  stable 
liquid  rocket  engines  without  excessive  trial  and  error.  One  such  empirical  result  is  the 
Hewitt  stability  correlation.  Figure  1  shows  the  correlation  for  rocket  engines  using 
LOX/RP-1  propellant  combinations  and  like-on-like  impinging-jet  injectors,  in  terms  of  the 
highest  sustainable  acoustic  frequency  and  the  injection  parameter  dJU-,,  where  is  the 
orifice  diameter  of  the  injector  and  U,  is  the  injection  velocity. 

Here,  we  are  concerned  with  identification  of  the  instability  mechanism  that  is 
responsible  for  the  emergence  of  the  Hewitt  correlation.  We  have  shown  that  such 
mechanisms  are  relevant  even  in  the  presence  of  turbulence  and  spatial  nonhomogeneities 
[1],  In  the  first  attempt  to  identify  the  main  acoustic  instability  mechanisms,  the  response 
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characteristics  of  various  combustion  processes  have  been  analyzed  to  calculate  their 
amplitude  and  phase  relationships  with  the  imposed  acoustic  waves.  Among  these 
processes,  the  dominant  amplification  effect  is  found  to  be  associated  with  near-extinction 
diffusion  flamelets  because  the  dominant  nonlinear  term  in  combustion  processes  arises 
from  the  finite-rate  chemical  reaction.  The  time-dependent  characteristics  of  near-extinction 
flamelets  have  been  analyzed  for  the  strained  diffusion  flame  model  by  employing 
activation-energy  asymptotics  to  predict  the  linear  response  of  the  burning  rate  [2,  3]. 
These  near-extinction  strained  diffusion  flames  are  expected  to  be  encountered  most  often 
between  impinging  spray  fans  of  fuel  and  oxidizer.  The  high  sensitivity  near  extinction  can 
be  seen  from  the  well-known  S-curve  behavior.  Fig.  2,  obtained  from  a  calculation 
performed  in  a  counterflow  diffusion  flame  [2],  shows  a  typical  result  of  the  steady- 
reaction- sheet  response  as  a  function  of  the  Damkohler  number  Da.  It  is  clear  from  Fig.  2 
that,  for  the  same  fractional  variation  of  Damkohler  number,  5  DalDa,  shown  by  vertical 
bars,  the  near-extinction  flame  experiences  a  much  larger  reaction-sheet  displacement,  and 
flame-temperature  variation,  as  marked  by  horizontal  bars,  also  is  much  larger.  In 
addition,  the  flame  responses  are  expected  to  diminish  as  the  acoustic  frequency  becomes 
greater  than  the  strain  rate,  since  accumulation  of  the  unsteady  effect  coming  from  the 
imposed  acoustic  oscillations  gives  rise  to  the  response  delay  and  attenuation  effects.  More 
detailed  results  on  the  flame  responses  with  acoustic  waves  are  discussed  in  our  published 
works  [2,  3,  4]. 

In  liquid-propellant  rocket  chambers,  the  rate  of  strain  decays  in  the  axial  direction. 
Consequently,  flamelets  near  extinction  may  be  expected  to  be  encountered  more  often  in 
the  near-injector  region  and  the  dominant  amplification  effect  is  considered  to  be  coming 
from  these  flamelets  near  the  injectors.  This  is  consistent  with  observations  in  engine  tests, 
indicating  that  stability  is  often  achieved  by  moving  the  flame  zone  away  from  the  injector 
assembly.  The  effects  of  finite-rate  chemistry  are  also  consistent  with  the  temperature¬ 
ramping  test,  a  stability-rating  method  in  which  the  fuel  is  gradually  heated  until  instability 
begins  to  disappear.  As  the  temperature  of  the  reactants  increases,  equilibrium  chemistry 
becomes  more  prevailing,  and  the  consequent  reduction  of  the  amplification  contributions 
results  in  stable  operation  of  the  rocket  engine. 

The  strain  rates  for  these  flamelets  can  be  approximated  as  U/d^,  which  is  the  strain 
rate  at  the  impact  point  of  the  two  jets  and  also  is  the  reciprocal  of  the  injection  parameter 
[5].  Since  the  strain  rate  in  impinging-jet  flows  tends  to  be  spatially  uniform,  the  value  of 
U represents  the  strain  rate  for  flamelets  in  a  relatively  wide  region  near  the  injector. 
However,  the  strain  rate  is  expected  to  fluctuate  about  U/d„  because  of  turbulence.  As  the 
strain  rate  U/d^  is  increased  by  increasing  [/,  or  by  decreasing  d^,  the  chance  to  encounter 
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near-cxtinction  flamelets  also  increases,  thereby  resulting  in  a  larger  amplification  rate. 
This  instability  tendency  is  consistent  with  the  Hewitt  correlation  in  that  instability  appears 
by  decreasing  the  injection  parameter,  i.e.  increasing  the  strain  rate.  In  addition,  the 
existence  of  a  high-frequency  cutoff  is  suggested  because  the  flame  response  rapidly 
decreases  with  increasing  acoustic  frequency.  A  substantial  increase  of  the  flame  response 
is  seen  in  our  previous  analysis  if  the  nondimensional  acoustic  frequency  co  is  smaller  than 
a  value  that  is  typically  less  than  unity.  Therefore,  the  cutoff  frequency  should  be  linearly 
proportional  to  the  strain  rate  UJd^,  which  is  again  consistent  with  the  Hewitt  correlation. 
Although  it  is  of  course  necessary  to  balance  amplification  against  attenuation,  and  many 
different  phenomena  arise,  it  is  interesting  that  this  one  particular  phenomenon  exhibits  all 
of  the  right  functional  dependences  [5]. 

4 .  Nonlinear  Acoustic  Response  of  LOX/GHj  Flames 

In  the  previous  year’s  summary,  some  of  the  results  on  the  linear  acoustic 
responses  of  subcritical  LOX/GH^  droplet  flames,  obtained  by  activation-energy 
asymptotics,  were  reported.  That  analysis  has  now  been  extended  to  include  the  effects  of 
detailed  H^-Oj  chemistry  and  nonlinear  acoustic  perturbations.  Since  the  pressure  is  well 
below  the  critical  pressure  of  oxygen,  the  characteristic  time  for  droplet  regression  is  much 
greater  than  the  characteristic  diffusion  time,  so  that  the  quasisteady  droplet  combustion 
model  is  still  valid,  and  acoustic  pressure  oscillations  can  be  supeiposed  on  the  quasisteady 
droplet  flames.  Also,  spatial  variations  of  the  acoustic  wave  are  not  considered  because  the 
characteristic  length  of  acoustic  waves  is  much  larger  than  that  of  droplet  flames. 

First,  the  steady-state  characteristics  of  the  droplet  flames  are  calculated  for  various 
pressures  by  using  an  inverse  numerical  method,  in  which  the  flame  structures  are 
calculated  as  a  function  of  the  maximum  temperature  instead  of  the  droplet  diameter.  Since 
the  flame  structure  is  a  single-valued  function  in  the  parameter  space  of  the  maximum 
temperature,  the  inverse  numerical  method  enables  us  to  calculate  the  flame  structure  even 
beyond  extinction  conditions,  where  the  conventional  flame  codes  fail  to  converge  because 
of  intrinsic  flame  instability.  Therefore,  the  inverse  numerical  method  provides  the  flame 
structures  near  extinction,  with  a  greater  accuracy,  where  the  flame  response  is  expected  to 
be  the  greatest.  Initial  results  of  the  steady  numerical  analysis  are  shown  in  Fig.  3,  where 
variations  of  the  maximum  flame  temperature  with  the  droplet  diameter  are  shown  for 
various  subcritical  pressures.  The  extinction  droplet  diameter  is  seen  to  increase  with 
decreasing  chamber  pressure,  and  extinction  can  be  achieved  for  a  droplet  flame  by 
sufficiently  reducing  the  chamber  pressure.  Therefore,  the  numerical  results  imply  that,  for 
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a  droplet  flame  with  a  given  diameter,  finite-rate  chemistry  is  more  important  for  low 
pressures,  while  equilibrium  chemistry  prevails  for  high  pressures. 

Nonlinear  response  to  the  imposed  acoustic  pressure  oscillation  P'  =  sin  (2ml x), 
with  the  acoustic  period  x  =10'^s,  is  considered  for  a  droplet  flame  of  2.5  /i/n-LOX- 
diameter  at  the  mean  pressure  P^  =  10  atm.  The  burning-rate  fluctuation  is  nondimension- 
alized  as 


sin 


^2m 

^  T 


\dt 


where  h  is  the  instantaneous  rate  of  heat  release  and  the  mean  rate  of  heat  release.  The 
variation  of  Q  with  the  normalized  acoustic  amplitude  PJP„  is  shown  in  Fig.  4.  The 
characteristics  of  the  heat-release  response  are  seen  to  be  extremely  nonlinear  in  that  the 
increase  of  Q  with  increasing  P^  becomes  greater  as  P^  increases  and  in  that  the  value  of  Q 
becomes  zero  beyond  a  critical  value  of  the  acoustic  pressure  amplitude  P^.  These  behaviors 
can  be  easily  explained  from  the  steady  flame  structure  shown  in  Fig.  3.  Since  higher 
flame  sensitivity  comes  from  the  effect  of  finite-rate  chemistry,  a  greater  nonlinear  heat- 
release  response  occurs  when  the  instantaneous  pressure  is  near  the  minimum.  As  the 
acoustic  amplitude  increases,  the  flame  becomes  closer  to  extinction  during  the  passage  of 
the  minimum  pressure.  Consequently,  the  increase  of  Q  is  faster  than^e  increase  of  P^.  If 
is  greater  than  a  critical  value,  which  is  about  4.9  atm  for  this  calculation,  then  the 
chemical  reaction  during  the  minimum  pressure  is  too  weak  to  support  the  flame,  so  that  the 
flame  becomes  extinguished  and  the  flame  response  vanishes. 

TTiis  t\pe  of  nonlinear  flame  response  is  viewed  as  a  possible  mechanism  to 
produce  the  threshold  phenomena  discussed  in  our  published  work  [6],  in  which  acoustic 
instabilities  can  be  triggered  only  by  acoustic  perturbations  with  amplitudes  greater  than  a 
critical  perturbation  amplitude.  To  show  this  possibility,  three  different  damping  response 
lines,  each  corresponding  to  strong,  weak  and  intermediate  damping,  are  drawn  as  dashed 
lines  in  Fig.  4.  Stationary  amplitudes  of  acoustic  waves  can  be  found  from  some  of  the 
intersections  of  the  amplification  and  damping  lines,  at  which  the  effects  of  damping  and 
amplification  can  balance  each  other.  First  considering  the  case  of  strong  damping,  the 
situation  marked  stable  in  Fig.  4,  the  only  intersection  occurs  at  P^  =  0,  so  that  any  acoustic 
perturbations  will  be  attenuated  and  the  system  is  stable.  For  the  case  of  weak  damping, 
there  are  intersections  at  P^  =  F,  and  at  P^  =  0,  and  only  the  intersection  at  P^  =  F,  is  stable. 
Therefore,  all  acoustic  disturbance  will  eventually  converge  to  F^  =  F^,  and  the  system  is 
always  unstable.  The  third  case  corresponds  to  intermediate  damping,  in  which  two  stable 
intersections  at  F^  =  P,  and  F^  =  0  are  separated  by  an  unstable  one  at  F^  =  F,.  If  an  initial 
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amplitude  of  disturbance  is  greater  than  P,,  the  acoustic  amplitude  will  converge  to  P,. 
Otherwise,  the  disturbances  will  be  attenuated.  Consequently,  the  system  is  metastable  and 
exhibits  threshold  behavior.  Our  analysis  has  demonstrated  that  such  a  metastable  acoustic 
system  can  exhibit  bimodal  distributions  for  the  amplitude  of  acoustic  waves.  It  is  also 
worthy  of  note  that,  though  the  nonlinear  response  shown  in  Fig.  4  is  calculated  from  a 
droplet-flame  model,  such  response  behaviors  are  expected  from  any  type  of  flame 
configurations  exhibiting  usual  high  sensitivity  near  extinction. 


5 .  Conclusions 

This  work,  in  total,  has  identified  a  particular  nonlinear  mechanism  involving  finite- 
rate  chemistry  and  potential  ignition  and  extinction  phenomena  that  can  lead  to  enhanced 
amplification  as  well  as  nonlinear  threshold  effects  in  acoustic  instability  of  liquid- 
propellant  rockets.  It  is  worthwhile  to  keep  this  new  phenomenon  in  mind  when  trouble¬ 
shooting  instability  problems  and  when  designing  motors  for  stable  operation.  Chemical 
rockets  do  involve  finite-rate  chemistry,  which  is  best  not  overlooked  in  addressing  engine 
designs. 

6 .  References 

1.  Kim,  J.S.,  “Effects  of  Turbulence  on  Linear  Acoustic  Instability:  Spatial 
Inhomogeneity,”  Liquid  Propellant  Rocket  Combustion  Instability  (V.  Yang  Ed.), 
Progress  in  Astronautics  and  Aeronautics,  169,  431-454,  AIAA.  Washington  DC, 
(1995). 

2.  Kim,  J.S.,  and  Williams,  F.A.,  “Contribution  of  Strained  Diffusion  Flames  to 
Acoustic  Pressure  Response,”  Combustion  and  Flame,  98,  279-299  (1994). 

3.  Im,  H.G.,  Law.  C.K.,  Kim.  J.S.,  and  Williams,  F.A.,  “Response  of  Counterflow 
Diffusion  Flames  to  (Oscillating  Strain  Rates,”  Combustion  and  Flame,  100,  21-30 
(1995). 

4.  Sohn,  C.H.  Chung.  S.H.,  Kim.  J.S.,  and  WiUiams,  F.A.,  “Acoustic  Response  of 
Droplet  Flames  to  Pressure  Oscillations,”  AIAA  Journal,  34,  1847-1854  (1996). 

5.  Kim,  J.S.,  and  Williams,  F.A.,  “Acoustic-Instability  Boundaries  in  Liquid-Propellant 
Rockets:  Theoretical  Explanation  of  Empirical  Correlation,”  Journal  of  Propulsion  and 
Power,  12,  621-624  (1996). 

6.  Clavin,  P.,  Kim,  J.S.,  and  Williams.  F.A.,  '‘Turbulence-Induced  Noise  Effects  on 
High-Frequency  Combustion  Instabilities,”  Combustion  Science  and  Technology,  96, 
61-84  (1994). 


5 


rfaqiiuocy  (Hoftz) 


4 


HHC  Csr- 


\  ^ 


0  HHC  H-t  ^envsiivs 
^  Lsf^C  !P3vli} 


3  Stable 
•  unstable 


r  H-1  §  0.73MN 

■‘S  O.SSMN^K 


5fablp  Grgion 


!  Cnstcbif  proicn 


HignesT  Sustainaole 
InstaDiliry  Fncusncy 


r  Orince  Diameter 
r.  Velocity  (seconds) 

::  suabiiiLV  correla:ion  for 
OX/PR-I  and  like-on-Iike 


Figure  2.  Steady  reac::oQ-5het 
name  tempera-ure  a5  a  inaction  of 
number,  sho^ag  the  high  seasitiv 
resDonse  near  e.xtinctioa. 


16  1^  1210  8 


Oroolet  ra(iius  [mm] 

ma/rm 

Figure  3.  Variaiioaof  the  maximum  temperature 
with  the  droplet  diameter  for  LOX/GK:  droplet 
hames  at  various  pressures. 

Figure  4.  Variation  of  the  nonlinear  heat-releas? 
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